ASYMPTOTIC BEHAVIOR OF QUANTUM WALKS ON THE LINE 
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Abstract. This paper gives various asymptotic formulae for the transition probability associated 
with discrete time quantum walks on the real line. The formulae depend heavily on the 'normalized' 
position of the walk. When the position is in the support of the weak-limit distribution obtained by 
Konno f |Kol| ). one observes, in addition to the limit distribution itself, an oscillating phenomenon 
in the leading term of the asymptotic formula. When the position lies outside of the support, one 
can establish an asymptotic formula of large deviation type. The rate function, which expresses 
the exponential decay rate, is explicitly given. Around the boundary of the support of the limit 
distribution (called the 'wall'), the asymptotic formula is described in terms of the Airy function. 



1. Introduction 

The notion of discrete time quantum random walks, often called quantum walks, was intro- 
duced by Aharonov-Davidovich-Zagury |ADZj in 1993 as a quantum counterpart of the classical 
1-dimensional random walks. Since then, quantum walks (and generalizations) have been inves- 
tigated in various contexts from both sides of mathematics and physics. See |Kej . |Ko2j for the 
historical backgrounds and the recent developments. In this paper, we shall give various asymptotic 
formulae for the transition probability associated with the 1-dimensional quantum walks. Our for- 
mulae explain the peculiar phenomena of the probability distribution on the whole interval [—1,1] 
(Figure [3Q) which is quite a bit different from the case of classical walks. 
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Figure 1. Probability disctribution and its weak-limit distribution 



The quantum walk we consider is concerned with a particle with spin 1/2 whose position is 
restricted to Z. To explain the set-up, we consider the Hilbert space 



f(Z,C 2 ) = {/ : Z 



C 2 
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with the inner product defined by 

(f,g) e := f(x), g(x) ), /,g G ^ 2 (Z,C 2 ), 

where || • || and (■, •) are the standard norm and inner product on C 2 , respectively. For x G Z and 
u G C 2 , define 5 X ® u G £ 2 (Z, C 2 ) by 

u (y = x), 



(8) w)(y) 



(y/x). 



Then, the set {5 X <g> ei,<5 x <g> e2}xez ; where {e\,e 2 } denotes the standard basis in C 2 , forms an 
orthonormal basis of £ 2 (Z, C 2 ). One has 

\\f(y)\\ 2 = \(f(y),e 1 )\ 2 + \(f(y),e 2 )\ 2 = \(f,5 y ®e 1 )e\ 2 + \(f 1 5 y ®e 2 ) £ \ 2 . 

Let A G SU(2), and write 

A=f^^, a,kC, |a| 2 + |6| 2 = l. (1.1) 



Decompose the matrix A aa3 



A quantum walk is described by a unitary operator £/ : £ 2 (Z, C 2 ) — )• ^ 2 (Z,C 2 ) defined by 

(Uf)(x) = P(f(x + l)) + Q(f(x-l)), f£l 2 (Z,C 2 ), xGZ. (1.3) 

For a vector 93 G C 2 with ||<^|| = 1, the quantity 

Pn(v;x,y) = \\U n (8 x ®(p)(y)f, (1.4) 

represents the transition probability that the particle with the initial state S x <8> <p is found at y 
after the n-step evolution. The unitarity of the operator U and the assumption that \\(p\\ = 1 tell 
that the sequence {p n (ip;x,y)} y£ z is a probability distribution on the lattice Z. As in the case of 
classical random walks on Z, the transition probability has the property 

p n (ip;x,y) =p n {tp;0,y - x), x,y£Z. (1.5) 

Thus it is enough to consider the case where x is the origin. Furthermore one can easily check that 
p n (ip;0,y) = when n + y is odd or \y\ > n. From now on, we set p n {(p;y) := p n (ip;0,y). Our 
aim is to give asymptotic formulae for p n ((p;y) as n tends to infinity. When a = or b = 0, the 
distribution p n (tp;x) can be easily computed. In fact, 

(1) if a = 0, we have p 2n (<p;x) = $x,o and p 2n+1 (<p;x) = \(fi\ 2 5 x _i + | ^2 1 2 ^,i , and 

(2) if b = 0, we have p n (f,x) = \ipi\ 2 5 x - n + I^P^.n, 

where <p = t ((pi,(p 2 ). Thus we only consider the case where ab 7^ 0. 

To state our main theorems, we define the function \a(<p) for (p G C 2 (||<^|| = 1) by 

\a{<p) = \<P2\ 2 ~ \<Pi\ 2 + T~^( abl Pi ( P2 + abtpiTp 2 ), ip = t {ipi,ip 2 ) £ C 2 . (1.6) 

M 

It is shown in [Kol] that the probability measure, 

dm nj(p = ^2p n (ip;y)5 y/n , (1.7) 



tSome of authors (as in |Ke] , |Kol| ) use the decomposition of the matrix A by the row vectors instead of the 
decomposition (|1.2[) . See below. 
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converges weakly to the one-dimensional distribution, 

X(-\a\,\a\)P(0^, p(£) = ^ i _p ) ^ s —^ , (1.8) 

where 8 y i n is the Dirac measure at y/n, and X(-|a|,|o|) i s the characteristic function of the open 
interval (— \a\, \a\). Roughly speaking, our asymptotic formulae for p n {ip;y) heavily depend on the 
normalized position y/n according as 

(1) y/n is inside the interval (— \a\, \a\), 

(2) y/n stays around the wall; say, y/n ~ ±|o|, or 

(3) y/n is outside the interval; say, \y/n\ > \a\. 

It is interesting to point out that our formulae described below resemble the Plancherel-Rotach 
formula for the Hermite functions (Theorem 8.22.9 in [Hi]), which tells us that the Hermite functions 
are 

(1) asymptotically oscillating when the position is in the classically allowed region, 

(2) approximated by the Airy function around the wall of the potential, and 

(3) exponentially decaying in the classically hidden region. 

In this view, 

(1) the interval (— \a\, |a|) is called the 'allowed' region, 

(2) ±|a| is called the 'wall', and 

(3) the outside of (— \a\, \a\) is called the 'hidden' region. 
The precise statement for the 'allowed' region is stated as follows. 

Theorem 1.1. Let a be a real number with < a < \a\. Then, for the integers y satisfying 

\Cn\ <a, i n = y/n, (1.9) 

we have 

p n (tp;y) = )J\ >H= L= [i + A ^)£ n + OSC n (£„) + 0(l/n)] (1.10) 

as n —7- oo uniformly in y satisfying (|1.9p . Here OSC n (£) is a function of the form 

A(£) cos(nfl(0) + B(£) sin(n0(£)). 

It is obvious that, for any £ G (— \a\, \a\), one can construct a sequence of integers {y n } satisfying 

y n = n£ + 0(l) (n^oo), (1.11) 

and the integers y in Theorem 11.11 can be replaced by such a sequence {y n }- Furthermore, if we 
allow {n} to be a sequence of positive integers, then we can construct a sequence {y n } satisfying 

y n = nti + 0{l/n) (n (1.12) 

by using the theory of continued fractions. For such a sequence satisfying (|1.12[) . all of £ n in (jl.lOp 
can be replaced by £. 

The factor l + (— l) n+Vn in the right hand side of (jl.lOp comes from the fact that Pui^', y) vanishes 
when n + y is odd. The functions A(^), B(^) and #(£) in OSC n (£) are computable (see Section [21 
in particular (|2.16p ). 

Observe that the distribution function p(^) defined in (jl.8p appears also in the asymptotic formula 
(jl.lOp with the factor 1/n. In fact, from Theorem 1 one can deduce the following corollary which 
is a special case of Theorem 1 in |Kol| . 
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Corollary 1.2. Let a, /3 be real numbers with —\a\ < a < f3 < \a\. Then we have 

r/3 



lim 

ySZ ; a<y/n</3 



MWV) = [ P(0d£, (1.13) 

^ , Jo 



where the function p(£) is defined in ([I 

The distribution p n (y?; y) mainly concentrates on the interval (— \a\n, \a\n) when n is large enough. 
However, p n {ip\ y) can be positive even when \y\ > \a\n. Hence it would be quite natural to consider 
its asymptotic behavior for {y n } with y n ~ ±n|a| or \y n \ > n\a\. 

Theorem 1.3. Suppose that a sequence of integers {y n } satisfies the following: 

y n = ±n\a\ + d n , d n = 0(n 1/3 ). (1.14) 

Then we have 

p n {<p; y n ) = (1 + (-l) n+ y-)a 2 n- 2 / 3 |Ai (ian' 1 / 3 ^) f (1 ± \a\\ A (<p)) + 0(l/n), (1.15) 
where a = (2 / 1 a. 1 1 6| 2 ) - 1 / 3 and Ai(x) is the Airy function. 

In the formula (|1.15|) . the variable in the Airy function is bounded, and hence its leading term 
is of order 0(ra~ 2 / 3 ). 

The asymptotic formula in the 'hidden' region takes the following form. 

Theorem 1.4. Let £ £ R satisfy \a\ < |£| < 1. Suppose that a sequence of integers {y n } satisfies 
(jl.lip . Then we have 

PniWVn) = { ] + { ~^ n ^^ e- nH ^\G(0 + O(l/n)), (1.16) 
7rn(l — £ )yc; — |a| 

where £ n = y n /n, G(£) is a smooth non-negative function in \a\ < |£| < 1, and i/ie function Hq(£), 
which is positive and convex in this region, is given by 

-OAna( \h\ 4- , /£2 _ U|2^\ 

(1.17) 



£T (0 = 2|e| log + V^ 2 - M 2 ) - 2 log (|6| + Ve^W 
+ (l-|£|)log(l-£ 2 ) -2|£|log|a|. 



This asymptotic formula tells that a large deviation property with the rate function Hq(£) holds 
in the 'hidden' region. Indeed, we have the following as a direct consequence of Theorem 11.41 

Corollary 1.5. Under the same assumption as in Theorem] 1.41 we hive 



lim - log p n (<p;y n ) = -Hq(£), (1.18) 

n— >oo n 

where n runs over positive integers such that n + y n is even. 

Here is a remark. As seen above, we have defined the unitary operator U by employing the 
decomposition (jl.2p in terms of the column vectors of A. Some authors use the decomposition 

a-* + s, *-(; J), ^(_° 5 1 

to define the unitary operator 

(Vf)(x):=R(f(x + l))+S(f(x-l)), fe£ 2 (Z,C 2 ), xG 
and the associated transition probability 

q n ^;x) = \\V n (6 ®^)(x)\\ 2 , VeC 2 , W = l. 
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It is easy to see that VA* = A*U on ^ 2 (Z, C 2 ), and from this we have 

q n (ip;x) = p n (Aip;x), x £ Z, ip£C 2 , \\tjj\\ = 1. 

Hence the use of different decompositions does not cause any significant differences in conclusion. 

We close Introduction by mentioning the strategy to prove Theorems ll.lt ES 11.41 and the orga- 
nization of the present paper. First, we express the transition probability p n (f,y) as the sum of 
the modulus square of the transition amplitudes; say, 

Pni^y) = \(U n (5 ®^), 5 y ® ei ) e \ 2 + \(U n (5 ®<p), 5 y <g>e 2 ) e \ 2 . (1.19) 

Thus it is enough to find asymptotic formula for the transition amplitudes ( U n (5o <g> </?), 5 y ® e% )t 
(i = 1,2). The starting point is the following integral formula 

1 

2iri 
-y 



(U n {5 ®<p),5y®il>} e = — I z-y- 1 (A(z) n v,i>)dz 



c 



2iri 



c 



(1.20) 

z~ y - 1 {A(uz) n <p,tj))dz, 



where uj = a/\a\ and tp £ C 2 with \\ip\\ = 1. The contour C is a simple closed path around the 
origin, and the matrix A{z) is defined by 

AM := P,- + Qz = ( a ±\ b j\ = C " !) ^I' °) . (1.21) 



6z y V— b aj V z 

The integral formula (|1.20p was also used in |GJS] to obtain the weak limit of the probability 
measure (jl.7p by considering the r-th moment of (jl.7|) . Our strategy is to analyze the integral 
(|1.20p by using a refined method of stationary phase due to Hormander (Theorem 7.7.5 in |Hoj ). 
Indeed, applying this method to the case where the contour C is the unit circle, one can prove 
Theorem ll.il (Section l2j). The proof of Corollary 11.21 involves a property of certain exponential sum 
determined by the asymptotic formula (jl.lOp (Lemma 12. 4p . Thus, we give the details of the proof 
of Corollary 11.21 in Section [2j When y/n stays around the wall, the phase function for this integral 
has degenerate critical points. But the third derivatives at the critical points do not vanish, and we 
can use the argument given in the proof of Theorem 7.7.18 in |Hoj to prove Theorem 11.31 (Section 
[3]). Finally, we change the contour C in the integral (11 .2Q[) to pick up a suitable critical point of 
the phase function and apply again the refined method of stationary phase. A detail of the proof 
of Theorem 11.41 is given in Section [H 

2. ASYMPTOTICS IN THE ALLOWED REGION 

This section gives a proof of Theorem ll.li To analyze the integral (|1.20p , we use a diagonalization 
of the matrix A(u)z). The characteristic equation of the matrix A(ujz) is given by 

A 2 - 2\a\(j){z)\ + 1 = 0, (2.1) 

where, for z £ C \ {0}, we set (f>(z) = (z + z~ l )/2. Thus the matrix A(poz) is diagonalizable except 
when the parameter z satisfies 4>{z) = ±l/|a|, or equivalently, when 

a= l±|6| 1±J6| 
\a\ \a\ 

Define the function A (z) by 

X(z) = \a\cp{z) + iy/1 - \a\ 2 (j)(z) 2 , (2.3) 
which is holomorphic on the domain 

D = C\{x €R; \x\ < (l-|6|)/|a| or \x\ > (1 + |6|)/|a|}. (2.4) 
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Then, for z G D, the eigenvalues of A(loz) are X(z) and X(z)~ 1 , where 



A(z)" 1 = \ a \(f){z) - iy/l - \a\ 2 (f>{z) 2 . (2.5) 

From ([23D, ([23D, we have X(-z) = -X(z)" 1 for z G D. Since (1 - |6|)/|a| < 1 < (1 + |6|)/|a|, the 
function X(z) is holomorphic around the unit circle. In the rest of this section, we take the unit 
circle S 1 for the contour C in the integral (|1.20p . Note that X(z) G S 1 for each z £ S 1 . Define the 
unit vector u(z) G C 2 for z G S 1 by 

«(*) = 4t f w P Z , . P(*) = V\b\ 2 + \Hz)~\a\z\ 2 . 
p{z) \X{z) - \a\zj 

For ■0 G C 2 with ||^|| = 1, define a vector G C 2 by 



The matrix [-0, J^] formed by two column vectors '0, is m SU{2). We may easily check 



A(uz)u(z) = X(z)u(z), A(ooz)Ju(z) = X{z)Ju{z). 

Thus any (p G C 2 is represented as ip = (tp,u(z))u(z) + (tp, Ju(z))Ju(z), and hence the integral 
formula (ll,20p becomes 

I _„_i . w 



(E^OSo® ®^>* = — / z ~ V ' lx ( z ) n M z ) dz + — [ . z-y- x X{z) g^{z)dz, (2.6) 
where we set 



27 ™ ./| 2 | = 1 27TZ J\ z \ = l 



f^{z) = (ip,u(z))(u(z),i)}, g^{z) = {ip,Ju{z))(Ju(z),il>). (2.7) 
Lemma 2.1. For any z G 5 1 , we have /?/,(— z) = g^{z). 

Proof. The vector u{— z) is an eigenvector of A{— ujz) = — A{ujz) with the eigenvalue A(— z). Note 
that, for z G S , we have A(— z) = — A(z) _1 = — X(z). Thus we see that A(uiz)u(— z) = X(z)u(—z), 
which means that u(— z) is an eigenvector of A{ujz) with the eigenvalue X(z). But, Ju(z) is also 
an eigenvector with the same eigenvalue X(z). Hence, there exists a constant c(z) G S 1 such that 
u{—z) = c(z)Ju(z), and we get 



U(-z) = ((p,c(z)Ju(z))(c(z)Ju(z),i(j) = c(z)c(z)(ip,Ju(z))(Ju(z),ij) = g^(z). 

□ 

On the unit circle, the function X(z) can be written as 

X(e u ) = \a\cost + - |a| 2 cos 2 t = e^ W , t G R, (2.8) 

where we set 

= Cos" 1 (|a| cost). (2.9) 

Note that < n(t) < vr. Define 

J(n,y) = J(^n,y) := -L / ^A^f/^) dz = ^- f e'^^ f^) dt. (2.10) 



2vri 7| Z | =1 2vr _ 

Then, from Lemma |2.1| it follows that 

( U n (5 ® A ® V> )/ = " -2/ (l + (-l) n+y ) n, I/). (2-11) 
Theorem 11.11 will be deduced easily from the following proposition. 
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Proposition 2.2. Define the function t(rj) on the open interval (— \a\, \a\) by 

t( V ) = Sin' 1 [ - p== ), N<|a|. (2.12) 
V |a| V 1 -»7 / 



Lei a be a real number with < a < \a\. Then, for the integers y satisfying (|1.9p . we have 

( U n (5 9 tp),5 y ® V >« = (1 + (-l)"+>^ L J® ^=== 

y 2vrn(l -OVl«l -Sn (2.13) 

X /ginM*($»))-Cnt(ln)]+Ti/ 4 ^( e i *(?n)-) _|_ e -«n[^(i(Cn))-€nt(Cn)]-^/ 4 ^(g-**(5n)^ -|_0(l/n)j 

uniformly in y satisfying (|1.9|) . 

Proof. For 77 with — |a| < 77 < |a|, we set 



Tr-tfa) (t?>0), 
-T-t(rj) (r ) <0). 

Note that, for 77 £ (— |a|,|a|), one has \t(r])\ < tt/2. Let e be a small positive number such that 
e < min{(7r-2t(a))/6, (\a\-a)/2, vr/8}. Take a function x(i) G C,f (-vr/4, vr/4) such that = 1 
for |i| < e and x(^) = for \t\ > 2e. We set Xn(t) = 1 — x(* ~~ *(£■«)) — ~~ *i(£n)), where, as in 
(jl.9p . we set £ n = y/n. Note that t (77) = ti(r/) if and only if 77 = ±|a|. The function x(t) is chosen 
so that Xn(t) = near t = i(£ n ) and i = ii(£ n ). The integral (|2.10p can be written in the form 

J(ip,n,y) = e~ iyt ^J (n) + {-\) n +v e iyt ^ n) Ji(n) + R(n), 

where Jo(n), J\{n) and R(n) are given by 

J (n) = — f e«M*+*««))-M/,( e i(W(en))) x ( t ) dtj 
2tt 

Ji(n) = ^f e- in ^^-*^» + ^^(e^*-*^)))x(t)dt, 

in|>(t)-£»t] f , 



2tt 

1 /•*Kn)+27T 



Note that, in the above, we have used 

A^t + tifa)) =7r-M(t-tfa)), /4e !(<+<l( " )) )=^(e !(H ' ,)) ), 

which is easily shown by using Lemma l2.1l and the definition of n{t). To use the method of stationary 
phase, we need the first and the second derivatives of fJ,(t), which are given by 

,. . lalsini „. . |o||6| 2 cost 

»'®= /t J -ttJ = ^' ^ * = 7TiM — 27w2- ^ 

y/l - \a\ 2 cos 2 t (1 - (a^cos^i)^ 

From (|2.14p . fJ,'(t) = 77 if and only if t = t{rf) or t = £1(77). Let cq > be the minimum of the 
function — rj\ on the complement of the set 

{(r?,i); |t-t(r/)| < e, |r?| < a + e} U {(77, t) ; \t-ti(rj)\ < e, \r)\ < a + e} 

in the rectangle [— \a\, \a\] x [— tt,tt]. Then we have \/jf(t) — £ n \ > cq near the support of Xn(t)- In 
the integral R(n), we extend Xn(t) as a function on K so as to be zero for £ > i(£ n ) + 2rr and for 
t < i(Cn)- Then it is still smooth depending on the parameter n and y. The support of the function 
Xn(t) on R, so obtained, is contained in the bounded interval [—t(a),t(a) + 2ir]. Note that, each 
higher derivative of Xn(t) is bounded uniformly in n and y, and hence we see R(n) = 0(n~°°) by 
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Theorem 7.7.1 in |Hoj . Since fjf(—t) = —/jf(t), it is clear that the functions /j,(t ± t(£)) =F £t have 
critical points only at t = in a neighborhood of the support of A direct computation gives 



V (t(v)) = V (-%)) = ^ , 

which shows that the critical point t = of the functions //(i ± *(£)) ~F & * s non-degenerate. The 
derivatives of functions f^(e l ^ t+t ^ n ^)x(t) and flty(e*(* - *« n ")x(£) are bounded uniformly in n, we 
can apply Theorem 7.7.5 in [Ho] to obtain 

J(ih,n,y) = — . 

V27rV'(*(^)) 

>ra[M*(&0)-£«*(f«)]+W4^( e »*(£n)) _|_ (_ l^n+J/g-m^^^n))-^*^)]-^/^ ^g-it(Cn) ) _|_ 0(1/ n) 

Note that, in the above, the term 0(1/ n) is uniform in y satisfying (|1.9p . Thus, multiplying 
the above by oj~ y (l + (—l) n+v ) and simplifying the terms involving (— l) n+y , we get the formula 
(12T3]) . □ 

Proof of Theorem[TJ\ Put /^(r?) = f ei (e u ^), g^ri) = g ei (e~ itM ) (i = 1,2), where {ei,e 2 } is the 
standard basis in C 2 . Define 

Taking the modulus square of (|2. 13)) for ij) = e\ and tp = e2 and summing up them, we get 



(l + (-i r+ ,)| b | 

( i + ( _i r+ ,)| 6 | 



Pn(¥\ y) = )l ' ~-f^== k„/l(Cn) + ^n5l(^n)| 2 + N„/ 2 (Cn) + ^n52(Cn)| 2 + 0(l/n) 



[|/l(^n)| 2 + |/ 2 (£n)| 2 + \gi(£n)\ 2 + |s 2 (£n)| 2 + OSC n (£ n ) + 0(l/n) 



7rn(l-e 2 )vW^ 
where the function OSC n (r/) is given by 

OSC„(7 ? ) = e 2m ^W^)-^)] + -/ 2 (/ 1 (r ? )^ + / 2 (r ? )^) 



+ e -^W*W)-^)]--/2( /l(?7 ) 5l(7?) + /2 ( 7? ) ff2(r?) ). 



I AMI 2 + 1/ 2 (??)| 2 + bi^)! 2 + \92(v)\ 2 = i + v[ M 2 - bil 2 + 7712(06^^2 + o6^ 2 ) 



Now, a simple, but a little bit long computation gives 

1 

\a\ 

which proves the asymptotic formula (jl.lOp . □ 
Remark: As stated in Introduction, the function OSC n (?7) can be represented in the form 

A(rf) cos(n%)) + B(rj) sin(n9(rf)), d(rj) = 2/j,(t(rj)) - 2rft(rj). (2.17) 

Proof of Corollary \1.SX Let a, (3 be real numbers such that — \a\ < a < f3 < \a\. For each positive 
integer n, we write 

{yo,Vl, ■ ■ ■ ,Vs n } = {y G Z; y + n e 2Z, na < y < n/3} (%-_ x < %■). 

Then, we have yj — yj-i = 2, j = 1, . . . , s n . Obviously s n = 0(n). For simplicity, we set 
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so that p(rj) = R(rf)(l + \a(<p)t]). By the formula ()1.10p . the sum in the left hand side of (|1 . 13)) is 
written as 

2 Sn 

M<P\ V) = -J2 R(VjM (1 + \A(<p)yj/n + OSC n ( yj /n) + 0(l/n)) . 



; a<y/n<P 



3=0 



The sum of the terms in 0(1 /n) is bounded because s n = 0(n). Hence, because of the term 2/n 
in front of the sum, the contribution from these terms vanishes when n tends to infinity. On the 
other hand, we find 

Z Yl + X A(f)Vj/n) = lim - ^ p(Vj/ n ) = / P(V) d-q. 

Tl—rOG (I Tl ?00 ffj I 

3=0 j=0 
In order to examine the sum ^ J2j=o R(yj/ n )OSC n (yj/n), we write 

J R(7 ? )OSC n (7?) = K(r])e in9 ^ +K^je- ind ^ l (2.19) 

where 9(rf) is defined in (12.170 and K(rj) is a smooth function on the interval (— \a\, \a\). Thus, 
Corollary 11.21 is a consequence of the following lemma. □ 

Lemma 2.3. Let K(rf) be a smooth function on [a,/3]. Then we have 

lim - K( yj /n)e ine{y i /n) = 0. 



n— >oo n 



(2.20) 



3=0 



To prove Lemma 12.3) we first show the following. 

Lemma 2.4. (1) Suppose that is not contained in [a, (3]. Then there exists a constant C > such 
that 

k 



3=0 



t in6(yj/n) 



< C 



for any k and n satisfying < k < s n . 

(2) Suppose that is contained in [a, 0\. Then there exists a constant C > such that 



£■ 



,inO(Vj/n) 



j=0 



< Cn l ' 2+5 



for any k and n satisfying < k < s n . 



Proof. First, we prove (1). We only consider the case a > 0, because the case /3 < is handled 
similarly. The first and the second derivatives of 9(rf) are given by 



?'(??) = -2*fa), 0"(v) = -2^(7?), 



(2.21) 



where t(rf) and R(rj) are defined in ()2. 12[) and (j2. 18[) . respectively. In particular, we have — tt < 
9'(rf) < on [a, 0\. The Taylor expansion gives us the equality 

e ine(y j+1 /n) = e in6(yj /n)+2i9' /n)+0(l/n) 
where 0(1/ n) is uniform with respect to j with < j < s n — 1. From this, it follows that 



Jn9(y j+1 /n) _ ind(yj/n) 



2i6'( yj /n) _ 1 



+ 0(l/n) (j = 0,..., Sn -l), 
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where we used the fact that the function f(jj) := l/(e 2ld '^ — 1) is bounded from above on [a,/3]. 
Using this formula and a partial summation, we have 

k fe-1 
^ e in%,/n) = ^ntffe/n) U{y ._ l/n) _ f{y ./ n)) + 

i=o j=l 

where the term O(l) is uniform in k with 1 < k < s n . Since f(rf) is smooth on [a,/3], we have 

\f{Vj-i/ n ) ~ f(Vj/ n )\ < C\Vj-i - Vj\/ n < C/n. 

Hence 



k 

ind(yj/n) 



3=0 



<Ck/n< Cs n /n < C, 



which shows (1). 

Next, let us prove (2). Note that we have 9'(0) = 0. The function 9(if) is written as 

9( v ) = 9(0) - r) 2 A( V ), A( V ) = -!{!- X)9"(X V ) dX. 

Jo 

Since 9"(r]) < 0, we have A(rf) > 0. We set f(rj) = rj^A(rf). Then, /'(0) = y/A(0) > 0, and 
hence one can choose e > such that f'(r]) > on (— e, e). Furthermore, since /(0) = 0, e can be 
chosen so that \2f(x)f'(x)\ < ir/8 on (— e,e). We assume first that the given interval [a,/3](9 0) is 
contained in (— e,e). We write 

k k—1 
^p e in9(y 3 /n) = ™ e »nfl(0) e -i(n 1/2 /fe/")) 2 (^+1 _ 

^ 2 ^ \ n n 

3=0 3=0 



We set Xj = f(yj/n). Since Xj + \ — Xj = 0(l/n), we have 

mi_y_i = r i (x j _ r i (x .) = 1 - x,) + o(i/n 2 ). 

n n /'(/ ^xj)) 

We further set tj = yfnxj. Then, tj+i — tj = 0(l/^/n), and 

— k-l 

n 



^yrrffo/n) = ^_ e m e (0 )jRn(A;) + ^ = ^ L ( tj / y^e^i (t j+1 - tj), (2.22) 

i=o i=o 

where = 1/ f'(f~ 1 (x)). Now, we claim that, for arbitrary 5 with 1/6 < 5 < 1/2, we have 

i? n (fc) = 0(n S ) (2.23) 

uniformly in k and n. For simplicity, we only consider the case where a < < f3, so that i Sn = 
y / nx Sn — >■ +oo and to = \/nxo — > — oo. Take a look at the sum 

S n (k):= ^ L(t j /y/n)e- i %(t j+1 -t j ). 

0<j<k-l ; \tj\>n s 

Note that when \tj\ < n s , we have = < n 5 " 1 / 2 . Since /(r/) is an odd function, we see 

\Vj/n\ < / _1 ( n5_1/2 ) = 0(n 8 ~ 1/2 ). Let l n be an integer such that yijn = (y + 2l n )/n = 0(l/n) 
as n — > oo. Then we have 

\(2j - 2l n )/n\ - |(y + 2/ n ,)/n| < | % -/n| < Cn 5 ^. 



ASYMPTOTIC BEHAVIOR OF QUANTUM WALKS ON THE LINE 11 

Therefore, we have \2j—2l n \ < Cn 5+l / 2 , which shows that the number of summands in R n (k) — S n (k) 
is of order 0(n &+l l 2 ). Since tj + \ — tj = 0(l/y / n), we obtain R n (k) = 0(n s ) + S n (k). Thus, it is 
enough to show that S n (k) is of order 0(n 5 ) uniformly in k. We shall only prove that the sum 

S n (k)= Yl Htj/^e^htj+i-tj) 

0<j<k-l;tj>n s 

is of order 0(n s ), because, for the sum in j with tj < —n s , the proof is carried out in the same 
way. We write 



e -^ 2 + i _ e -«? = e ~ it2 j ( e -*y _ i) j a . = 2tj(t j+l - tj) + {tj+i - tj) 2 = t 2 j+1 - t 



Since tj = ^frixj = \/nf(yj/n), we see 

tj(t j+1 - tj) = nxj(x j+ i - Xj) = 2f'(yj/n)xj + 0(l/n) = 2f(y j /n)f'(y j /ri) + 0{l/n). 

Thus a,j = 4:f(yj/n)f'(yj/n) + 0(l/n). Let j n be the smallest number such that tj n > n s . Since 
< 2f(x)f'(x) < it/8 on (0,e), we have < aj < 7r/4 for n large enough, and hence e m i ^ 1 for 
any j n < j < s n . We shall check that 

^-^ = 0(n^), j = j n ,...,s n . (2.24) 

Indeed, we have 

\ e ~ ia 3 — 1| > sinaj = ajAj, Aj = / cos(ajA) dX. 

Jo 

Since < aj < it/4, there exists a constant c > such that Aj > c. Thus, for j n < j < s n , 



e -ia 3 _ i 

which shows f|2.24j) . Substituting now 



1 1 A 

< -7 r < — n , 

c[tj + i + tj) 2c 



-it 2 e 3+i - e j 



e~ w i = : (2.25) 

for S n (k) and using a partial summation, we get 

fc-i 

S n (k) = ( D J-i ~ D j> + D k . ie ~ u l + D jn e' u l , 

j=jn + l 

where we set 

Dj = t J^L L (t j /y^). (2.26) 

In view of (|2.24|) . we see 

fe-i 

Sn(k)= Y {Dj-i-D 3 )e~ it2 i +0(n~ s ). (2.27) 

j=jn + l 

Since L(x) is a smooth function and since tj + \ — tj = 0(1/ \/n), we have 

Dj _ x - Dj = L(t 3 /VTi) (^r^Y " ^7^) + 0(n~^). (2.28) 
To estimate the right hand side of (|2.28p . we write 
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Since 

t j+1 = Sfif(y j+1 /n) = tj + 2f'( yj /n)/^ + 0(l/n 3/2 ), 
tj- x = tj - 2f'( yj /n)/^ + 0(l/n 3 / 2 ), 

we have 2tj — tj + x — tj-x = (9(l/n 3 / 2 ). ()2.30p also shows that tj + x — tj > c/y/n with a constant 
c > 0. Thus we have l/(e _4Clj — 1) = C^n 1 / 2-5 ), and hence the first term in the right hand side of 
(ICTjl is of order 0{n~ l ~ s ). Note 

0(n~ 2S ). (2.31) 



g-lOj-l _ X e ~ ta ] - 1 

Indeed we write 

1 1 e - ia i-i 



^ e l(Oj_l-Oj) _ _ 



we 



g-ifij-i _ x e - io j — 1 (e~ ia i- 1 — l){e~ ia i — 1) 
Setting h{x) = f(x) 2 , we have 

t 2 j+1 = nh(y J+1 /n) = t 2 + 2h'( yj /n) + 0(l/n), t 2 ^ = t 2 - 2h'( yj /n) + 0(l/n). 

From this we have aj - a 3 -x = t 2 j+1 + t 2 _ x - 2t 2 = 0(l/n). Since l/(e - ^ - 1) = C^n 1 / 2 " 5 ), 
obtain (|23T|) also shows that ^ - = O^" 1 / 2 " 25 ). Therefore, by (l237jl . we have 

S n {k) = 0(n 1//2_2<5 ). But since 5 > 1/6, we obtain S n (k) = 0{n 5 ) as required. 

Next, we shall prove (2) for general [a, f$\ containing 0. As before, we may suppose that a < < 
p. Take c > such that c < min{|a|,e, f3}. We write 

k 

y^ginflfe/n) _ e in9 (Vj/ n ) _|_ ^ e in6{Vj/n) _|_ e in9{ yj /n) 

3=0 0<j<k;yj/n€[a,—c] 0<j<k;yj/ne[c,/3] 0<j<k;-c<yj/n<c 

By (1), the first and the second sums are of order O(l). Since each term of the sum is 0(1), we 
have 

e ind(yj/n) _ e ine (y 3 / n ) _|_ O(l). 

0<i<& I —c<yj /n<c 0<j<k ; —c<yj /n<c 

By the fact we have just proved, the sum in the right hand side above is of order 0(n 1//2+5 ), and 
hence we obtain the desired estimate. □ 

Proof of Lemma \2.cH We set 

k 

T fc (n)=r fc (n;[a,/3]):=^e m ^/"). 

j=0 

If [a, 0\, then Tfc(ra) = 0(1) uniformly in k. which is, by Lemma 12.41 of order 0(n 5 ) for any 5 
with 1/6 < 5 < 1/2 uniformly in k = 0, 1, . . . , s n . A partial summation leads us to 

- £ (y>)e in ^>) = - J] 7» - K(y j+1 /n)) + 0(n^+ 5 ) 

71 j=0 71 j=l 

Since ^(77) is smooth, we have K(y 3 /n) — K(yj + x/n) = <9(l/n). Thus, Lemma [23] shows that the 
first term in the above is of order 0(s n n s ~ 2 ) = 0(n~ 1+s ). Therefore, we conclude (|2.20p □ 



(3.1) 
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3. ASYMPTOTICS AROUND THE WALL 

In this section we give a proof of Theorem II .31 As in the previous section, we take the unit circle 
S 1 for a contour C in fll.20f) . However, the phase functions //(i)^ H^> which appears in the integral 
(|2.10p . have degenerate critical points, and thus we cannot use Theorem 7.7.5 in [Ho] . Instead one 
could use Theorem 7.7.18 in [Ho] directly to obtain the following proposition. 

Proposition 3.1. Suppose that a sequence of integers {y n } satisfies (|1.14p Then we have 

(U n (5 ®<p), S yn <S)^) e 

= w-v»j*("*vn)/2(i + (-l)"+^)/^(±i)an- 1 / 3 Ai(±an^ 1 / 3 (i n ) + 0(n~ 2 / 3 ). 

where a = (2/|a||6| 2 ) 1//3 , and Ai(x) is the Airy function defined by 

Ai (x) = — / c'^ 3 +^d£. 
2vr J_ O0 

Proof. For simplicity, we consider only the case where y n = \a\n + d n with the assumption d n = 
0(nV 3 ). -^/g neec i t analyze J(ip,n,y n ) defined in (|2.10j) . For this sake, take a function p(t) £ 
Cq°(— 7r/4, 7r/4) satisfying p(t) = 1 near t = 0. Set pi(t) = 1 — p(t), and write 

Jty,n,y n ) = Ji(n) + J 2 (n), J x (n) = — / e m ^ M p{t) U{ie~ lt ) dt, 

27r J-tt/2 

where the integral J2{n) is defined by replacing p(t) with pi(t) in the definition of J\(n), and the 
function <p is given by 

(f)(t) = p(-t + n/2) + \a\t. 
The derivatives of (j> up to the third order are given by 
, \a\cost „ |a||6| 2 sint ,„ \a\ \b\ 2 cos t{l + 2|a| 2 sin 2 1) 

( P[t) = -T l i !9 ■ 2,M/9 +l°l' ^ (*) = 7i 1 19 ■ ^ (*) ~ 



(1- laPsin 2 *) 1 ^ 1 " rw (l-| a | 2 sin 2 t) 3 /2' r w (1 - |a| 2 sin 2 1) 5 /2 

Thus, on the interval [— tt/2, 37t/2], <j>'(t) = if and only if t = 0, and we have 

0(0) = tt/2, 0'(O) = 0"(O) = 0, 0"'(O) = |a||&| 2 > 0. 

Since <j>'(t) / on the support of pi(t), a standard argument involving an integration by parts 
shows J2(n) = 0(n~°°). One may also write 



4>(t) = 7r/2 + t 3 A(t), A{t) = ±J(l-x 



2 4>"'(xt) dx. 



Note that A(0) = |a||6| 2 /6 > 0. Set B(t) := (3A(t)) 1 / 3 t on a small neighborhood of t = 0. Then 
we have 

<j ) (t) = Ti/2 + B{tf/? > , 5(0) = 0, ^'(0) = a" 1 = (|a||6| 2 /2) 1 / 3 > 0. (3.2) 

Choosing p{t) whose support is small enough, we may suppose that B'(t) > near the support of 
p(t). A change of variable leads us to 

*(„) = ^1 /e<»<~ -,„/»,„„ ds , u(s) = Wpll . ,3.3, 

Taking /o(i) whose support is small enough if necessary, we may suppose that u(s) is a compactly 
supported smooth function. One can now apply Theorem 7.7.18 in |Hoj directly to the integral 
J\[n) to finish the proof. For the sake of completeness, we shall give a detail. Define a function 

(f(s,r) by 

tp(s,r) = s 3 /3 + rB~ 1 (s). (3.4) 
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Since 

<p(o,o) = (d s <p)(o,o) = (d s V>(o,o) = o, (a s V)(o,o) = 2, 

we can apply Theorem 7.5.13 in |Hoj which tells that there exists smooth real- valued functions 
T(s,r), a(r), b(r) defined on a neighborhood of (s,r) = (0,0) satisfying 

a(0) = 6(0) = T(0, 0) = 0, (d s T)(0, 0) > 0, 

1 , (3-5) 

ip(s, r) = -T(s, rY + a(r)T(s, r) + b(r). 

By (13. 4p . the functions T, a, 6 have the properties 

T(s, 0) = s, a'(0) = a, 6'(0) = 0. (3.6) 

Indeed, setting r = in the second line of (|3.5p . we have T(s,0) = s. Computing the derivative 
(d r ip)(s, 0), we also get 

= 6'(0) + a'(0)s + s 2 (d r T)(s, 0). 

Since -B _1 (0) = 0, (-B _1 )'(0) = a, we have (|3.6p . By the implicit function theorem, one can find 
a smooth function s(t,r), defined on a neighborhood of (t,r) = (0,0) such that T(s(t,r),r) = t. 
Since T(a,0) = s, we have s(0,0) = 0, (<9 t s)(0,0) = 1, (<9 t 2 ,s)(0, 0) = 0. A change of variables in the 
integral (j3.3|) gives 

«7r(n— j/n)/2 /■ 

Jl(n) = 2^ / em(i3/3+a(/3n) * +6(/3n)) ^(^/ 3 n)^ ) = r))(5ta)(t,r), (3.7) 

where we set j3 n = d n /n which is of order 0(n _2//3 ). Choosing p(t) whose support is small enough, 
one can assume that the function v(t,r) is compactly supported around (t, r) = (0,0). Thus we 
can choose a function %(t) £ Cq°(R) such that < x(t) < 1 an d x(0 = 1 when v(t,r) ^ 0, and 
hence v(t,r) = x(t)v(t,r). The Malgrange preparation theorem (Theorem 7.5.6 in |Hoj ) tells that 
there are smooth functions q(t,r), c(r), d(r) satisfying 

v(t, r) = q(t, r)(t 2 + a(r)) + c(r)t + d(r). (3.8) 

Note that v(0, 0) = u(0), (dtv)(0,0) = u'(0) where the function u(s) is defined in (13. 3|) . Thus we 
have 

d(0) = u(0), c(0) = u'(0). (3.9) 
Inserting (]3.8h into (|3.7p and integrating by parts, we obtain 



„iir(n—y n )/2 



+ c(A0 / e^/^^txC*)^- — / e in (* 3 / 3+a ^)*)^(x(t)g(t J( 5n))^ 



The last term in ()3.10|) is of order 0(l/n). For the first and the second integrals, a change of 
variables shows 

e m(*V3+a(/3„)t) x ( t ) ^ 

2vr ? i- 1 /3Ai(n 2 / 3 a(/3 n )) - n" 1 / 3 />* 3 /3+n 2 / 3 a(^)*) (1 _ x { n - l 'h)) dt, 

J (3.11) 

e in(t»/3+«(A,)t) tx(t ) df 



-27rm- 2 / 3 Ai'(n 2 / 3 a(/3 n )) - n~ 2 / 3 | /^MMt)^ _ x („-i/3 t) ) 



dt. 
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By the assumption j3 n = d n /n = 0(n _2//3 ) and the properties (I3.5p . (I3.6p . we have 

Ai(n 2 / 3 a(/3 n )) = Ai{n 2/3 ap n ) + 0(n 2 / 3 /? 2 ) = A\{an~ l ' z d n ) + 0(n" 2 / 3 ). (3.12) 

By the same reason, the term n _2//3 Ai'(ri 2 / 3 a(/3 n )) in the second line of (|3.1ip is of order 0(n~ 2 / 3 ). 
To estimate each integral in the right hand side of (I3.11D . we need the following lemma. 

Lemma 3.2. Let a n be a sequence of real numbers such that a n — >■ as n — > oo. Let < 5 < c, 
e > be constants satisfying 5 < 1, 5 + e < c — e, and x( s ) be a smooth function on R such that 
x(s) = 1 for \s\ < 5 + e and x(s) = for \s\ > c — e. Then the integrals 

e in(s 3 /3+a n s)^ _ ^ J ^/S+Ons) g ^ _ dg (3^3) 

are of order 0{n~°°) as n — >■ 00. 

Assuming Lemma 13.21 one sees that the integrals in the right hand side of (13. lip is of order 
0(n-°°). Using d{p n ) = u(0) + 0(n" 2 / 3 ) and nb((3 n ) = 0(np 2 n ) = C^n" 1 / 3 ), we get 

p iTr(n—y„)/2 . . 

j^n) = e inb ^ ( 2vrn- 1 / 3 n(0)Ai(Qn- 1 /3 ( i n ) + 0(n" 2 / 3 ) ) 

= e iw(n - y " )/2 n- 1 / 3 u(0)Ai(an- 1 / 3 d n ) + 0(n" 2 / 3 ). 
By the definition (|3.3p of u(s), we have u(0) = af^(i), and hence we conclude the proof. □ 

Proof of Lemma \3.2\ For the integrals f)3. 13j) , the integrals on the negative real line have the same 
form as the integrals on the positive real line except the sign of n. Thus it is enough to prove that 
the integrals 

roo poo 

#i(n):= / e in ^ s Xl- X {s))ds, R 2 (n) :=/ j n ^s(l - X (s)) ds 
Jo Jo 

are of order 0(n -00 ), where we set VVi(s) = s 3 /3 + a n s. We give a detail only for the integral R%{ri) 
since the arguments for R\{n) is the same. We may assume \a n \ < 5 2 /2 for every n. Take R > c 
and write R2{n) = S(n, R) + T(n, R), where 

pR poo 

S(n,R)= / e in ^^s{l - x (s))ds, T(n,R)= / e in ^sds. 
Js Jr 

since ijj'(s) = s 2 + a n / for 5 < s, we can use the differential operator tfi'(s) 1 {d/ds) to perform 
an integration by parts, and get 

in(R z + a n ) m Js s z + a n inj s (s z + a n ) z 

where we have used that X (R) = and X (S) = 1. By the assumptions of X (t), the function 
d s (l — x( s )) = ~x'( s ) nas a compact support on the interval (5,R), and hence the integral in the 
second term is of order 0(n~°°). Repeating this procedure, we have, for any N > 0, an expression 
of the form 

S(n R) = V ^^ (fl> +^-[* e^Wfl - X(s)) ds + 0(n-) 

where qk{s) is a polynomial of degree at most k independent of R. Note that the estimate 0(n~°°) 
does not depend on R (but depends on N). From this expression, it is easy to show that 

S(n,n N - 1 ) = 0(rT N ). 
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Next, let us consider the integral T(n, R). We set /(C) = (e in ^ 3 / 3+an ^ , ( S C. We take L > and 
77 > such that R < L and < rj < c/2. Integrating /(C) on the rectangle determined by R, L, 
L + ir], R + irj and letting L — > 00, we have 

pOO pTJ 

T(n,R)= f(s + ir])ds + i f(R + it)dt. (3.14) 
Jr Jo 

It is clear that 

\f(R + it)\< CRe~ n ( R2t - t3 / 3+a ^ < CRe~ nRH/2 , < t < 77, 

where C > is a constant independent of n, R, and R > is chosen so that R 2 — 5 2 > 2rj 2 . Thus 
the integral in the second term in (|3.14p is of order 0(1/R). The first term in (|3.14p is estimated 
as 

poo 

\f(s + ir])\ds<C / se - n n{s 2 -r?/Z+a n ) ds < Cn -l e -nr,(R?- V y3+a n )_ 

Jr. 

Thus, if we take R = n k with k > 1, then the first term in (|3.14p is of order 0(n~ 1 e~ vn2k+1 ^ 2 ), and 
hence T(n, n^" 1 ) = 0(n~ Ar+1 ). Therefore, we get 

R 2 (n) = S(n, n N ) + T(n, n*) = 0(n~ N ) 

for any A^ > 0, which gives the assertion. □ 

Proof of Theorem \1.3l Using (|2.7p , we can easily check 
and hence we have 

|2 , I, /,-m2 1 , l°l A |2 1 |2 , 06^(^2 + O Vi^ 2 \ 1 



i/«i(±»-)r + i/ e2 (±*)f = 2 ± y (j^ - bir + , a|2 ^^ j = §a ± maaH)- 

Taking the modulus square of (|3.ip for ip = e\,e2 and summing up them, we conclude (|1.15p . □ 

4. ASYMPTOTIS IN THE HIDDEN REGION 

The purpose of this section is to prove Theorem 1 1.41 As in the previous sections, we first handle 
the transition amplitude. 

Proposition 4.1. Let £ satisfy \a\ < |£| < 1, and suppose that a sequence of integers {y n } satisfies 
(jl.ll|) . Then we have 

e -nH Q (tn)/2 



-!Jn 



\/2im 

(4.1) 

b\ 



a. 



^ n - M/2 (F^) + 0(l/n)), 



where £ n = y n / n o-nd -^(C) is a smooth function in \a\ < |£| < 1 which is given in (|4,36p . (|4.37p . 
and Hq{^) is a convex positive-valued function defined in (|1.17p . 

In the previous sections, we have taken the unit circle for a contour C in the formula (|1.20p . 
The phase function, however, has no critical points when the parameter £ is outside the interval 
(— \a\, \a\). Therefore, it would be reasonable to change the contour C to pick up a suitable critical 
point. Here a 'phase function' means essentially the logarithm of the eigenvalues of the matrix 
A(uz). 
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Before proceeding to the details, let us give a rough description of the phase function and its 
critical points. Recall that, when z £ C \ does not coincide with the values in (I2.2D . the matrix 
A(ujz) is diagonalizable. Hence, by diagonalizing A(uiz), (|l,20p is expressed as a linear combination 
of integrals of the form 

" z-y- 1 X(z) n f(z)dz, 
c 

where f(z) is holomorphic on a domain containing C, and X(z) is an eigenvalue of A{ujz) (we 
assume here that \{z) is also holomorphic). By expressing the path C by a function z = z(t) with 
z'(t) 7^ 0, the integral above is written as 

e n ^f(z(t))^ dt, $(t) = logA(z(t)) - 7/logs(t), 

where we set r/ = y/n. The critical point c of the 'phase function' satisfies the equation, 

z W)~ r, = {) {z = z(c)) - (4 - 2) 

Our strategy is to change the path C so that it passes through a point satisfying the equation f|4.2j) 
and is to use Theorem 7.7.5 in |Hoj . To make it possible, we need to cut the complex plane so as 
to be able to take the eigenvalue \(z) depending holomorphically on z. 
In what follows, we divide the discussion into the following three cases: 

(A) - 1 < |£| < 1, (B) \a\ < |e| < ; 1 (C) |e| = ; 1 (4.3) 

In the rest of this section, the following two functions play important roles: 



m = \n+jrw t rfa) = (W<H<1) (44) 

yjl-rf \a\s/l-r] z 



Note that one has 



Djrf)-^- l b l ~ V 7 ? 2 ~ l Q l 2 r(r?) -i _ I 6 ! 7 ?" V^-I^l 2 



a/1 - ?/ 2 l"lvl - ^ 

The function Hq{t]) defined in (|1.17j) can be expressed as 

H Q (r])/2 = 77 log \r(rj)\ - logD(rj) = \r]\ log r ( 1 77 1 ) - log D(rj). (4.5) 

4.1. Proof of Proposition 14.11 for the region (A). In this subsection, we assume that the 
parameter £ in the statement of Proposition 14.11 is in the region (A) . 



4.1.1. Some lemmas. Let Vm(C) De a function in £ £ C \ [—1, 1] defined by 

^a(C) = C + (C-i)< 



'C + i 



. - r 

where the branch of the square root is chosen so that it is positive on the positive real axis. The 
function (j){z) = (z + z~ l )/2 maps the domain {\z\ > 1} conformally and injectively onto C\ [— 1, 1], 
and the function Vm(C) is its inverse. Thus, if we set 

\a(z) := ipA(\a>\(/>(z)), z £ D A , 

then the eigenvalues of A(luz) are given by \a{z) and \a{z) _1 , where the domain Da is defined by 

D A := C \ (S 1 U {0} U [-(1 + \b\)/\a\, -(1 - |6|)/|a|] U [(1 - |6|)/|a|, (1 + |6|)/|a|]) . 
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The function \a(z) is obviously holomorphic on the domain Da- Define the vectors ua(z),va(z) E 
C 2 by 

(ab \ / ab \ 

, n R , -i , v A {z)=[. F .J. (4.6) 
Aa(z) - |a|z V \Aa(z) - |a|z V 

One can easily check that the vector ua{z) (resp. va(z)) is an eigenvector of A(ujz) with the 
eigenvalue Aa(^) (resp. Aa(-z) -1 )- For any cp = (971,(^2) G C 2 , we set 

|a|(|a|^ -1 - Aa(2;) _1 )^i + a6zy> 2 |a|(|a|z _1 - A j4 (^))v?i + ^ 

z = h t\ t \ \ 1 \-\\ ' <P2A(z) = 4.7) 

abz{\A{z) - \a{z) l ) abz{\ A {z) - \a{z) l ) 

so that 

= 9?ia( 2 ) u a(^) + ^2a(z)va(z), ^ e -Da- (4.8) 

For ^ £ C 2 , we set 

fAi/>(z) = Via{z)(u a {z),iIj), Qa^z) = ip2a(va(z),^}. (4.9) 
Then the integral formula (|1.20j) is written as 

( U n (5 Q ®<p),8 v ®il>)t = uj- y {I A i{n) + I A 2(n)), 



lAi(n) = ^- [ z~y- l \ A {z) n f AiJ (z)dz, I A2 (n) = ^- [ z'^ 1 X A (z)- n 'g A ^(z) dz. 
2tti J c 2-ki J c 



(4.10) 



Lemma 4.2. Let rj € R with \a\ < \r]\ < 1. Then the solutions of the equation (|4.2p with A(z) 
Xa(z), z S Da, are given by z = ±ir(r]), where r(rj) is defined in (|4.4p . 

Proo/. Since </>a(C) = J ^§V(C - 1)/(C + 1), w e have 



A^z) |q|/i(z) / |q|^)-l _ z- z 1 

^Aa(^) " H0(z)-lVH^) + l' M j " 2 ■ l4 - iij 

From this and ^(z) 2 — fJ-(z) 2 = 1, we see that the solution z 6 Da of the equation (|4,2p must satisfy 
(p(z) 2 = — (rf — |a| 2 )/|a| 2 (l — rf). Hence z is of the form z = iy with a real number y. Since 

<fi(iy) = ifJ-(y), we see fj,(y) = ±-\A? 2 ~~ l a l 2 /l a l \A ~ V 2 - Since 

ir + 1 \-Mf ( r< °)' 

we have 4>(y) = ±|6|?7/|a|\/l — rj 2 . Our claim is a consequence of the relation y = fj,(y) + 4>(y)- □ 

The following lemma is used to find a suitable contour C in (|4.1U|) . 
Lemma 4.3. Let r > (1 + |6|)/|o| or < r < (1 - |6|)/|a|. T/ien we have 

|Aa(±0| < |AA(re l *)| < |A A (±ir)| (-tt < * < tt). 

The equality in the first inequality holds if and only if t = 0, ±ir, and the equality in the second 
holds if and only if t = ±tt/2. 

Proof Let r > (1 + |6|)/|a|, and let C r be the circle of radius r centered at the origin. The image 
of C r under the map 4>(z) = (z + z _1 )/2 is the following ellipse 

4>(C r ) : 4>(re lt ) = <p(r) cost + ifi(r) sint (-tt <t <tt), 

and thus the curve \a\4>(C r ) is also an ellipse, where, as before, we set fi(z) = (z — z~ l )/2. Take 
two positive numbers x, y such that 

4>(x) = \a\4>(r), fi(y) = \a\fj,(r), x,y > 1. 
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The numbers x, y satisfying the above are uniquely determined. The ellipse \a\4>(C r ) lies outside 
4>(C X ) and inside 4>(C y ). Furthermore \a\4>(C r ) is tangent to 4>(C X ) at the points ±|a|0(r) = ±<f>(x), 
and tangent to 4>(C y ) at the points ±i\a\/i(r) = ±i/x(y). Hence the curve A J 4(C r ) = ijjA(\a\4>(C r )) 
lies between the two circles C X} C y . Since the map Vm is conformal, the curve \a{C t ) is tangent 
to C x at ±x from the outside of C x and is tangent to C y at ±iy from the inside of C y . Thus we 
have x < |A^(re**)| < y. Since 4>(iy) = i/J-(y), we see 

±iy = VuO(±*y)) = ipA(±iv(y)) = ^A(±i\a\n(r)) = ip A (\a\(j)(±ir)) = Aa(±«>), 

and similarly ±x = A^(±r). This proves the desired inequalities. The conditions for the equalities 
are a consequence of the fact that the curve A J 4(C r ) is tangent to C x only at ±x and to C y only at 
±iy. Next, let us consider the case where < r < (1 — |6|)/|o|. Since |a|/(l — \b\) = (1 + |6|)/|a|, we 
have (1 + |£>|)/|a| < r _1 . Thus, from what we have just proved for the case where r > (1 + |6|)/|a|, 
we see 

lA^ir" 1 )! < IXAir-'e^l < \X A (±t r - x )\. 

Noting 4>{z~ x ) = 4>{z), we see \a{z~ 1 ) = Xa{z), and hence replacing t by —t in the above inequality, 
we obtain the desired inequalities together with the conditions for the equalities. □ 



4.1.2. Proof of (|4.ip for the case £ > in the region (A). The condition 1/a/1 + \b\ 2 < £ < 1 is 
equivalent to the inequality (1 + |6|)/|a| < r(£). Set £ n = y n /n as in the statement of Proposition 
14.11 and take the circle C n = CV(£ n ) of radius r(£ n ) for the contour C in (I4.10|) . Since £ n , — > £ as 
n —> oo, we have r(£ n ) > (1 + |b|)/|a| for every sufficiently large n. The integral lAi(n) with y = y n 
defined in (|4.10p can be written in the form 

e -n[£„logr-(£ n )-logD(£ n )] 



lAi(n) 



2vr 



-iy n t 



e -nff Q (e„)/2 


f'7T 

/ e -iynt 


[A A (r(£ny*)l 


2vr 


— 7T 





n 

fAMZnY^dt 



fAip(r(tn)e n ) dt 



where -D(£) is defined in (|4.4|) and we have used (|4.5p . A direct computation shows 

A^(±ir(r/)) = ±iD(r]), \a\ < \rj\ < 1. 



(4.12) 



By Lemma 14.31 for each 77 with 



< 77 < 1, we have the inequality \\A(r(r])e lt ) / D(r])\ < 1, 



and the equality holds if and only if t = ±7r/2 on the interval [— vr,7r]. Take a function £ 
C^°(-7r/4,vr/4) with x(*) = 1 near t = 0. Setting xi(i) = 1 - x(t - 7r/2) - %(t + tt/2), the integral 
Iai(") can be written as 

J A1 (n) = e -^(«")/ 2 (A + (n) + A_(n) + R(n)), 

where we set 

=™/ 2 r . . rA^(±ir(e„)e li ) 



A±{n) 



R{n) 



2ir 
1 

" 2vr 



- ilfe V^(±ir(£ n ) e ")x(t) 



" l?/ ' lt /A^(r(en)e J ')xi(t) 



AA(r(^„)e 4 *) 



dt, 



dt. 



We can choose a positive constant < c < 1 such that 

AA(r(^)e 1 *) 



< c 



for every n and every t near the support of Xi(t). Setting R = — log c > 0, we see = 0(e nR ). 

Next, take a look at the integrals A±(n). By (|4.12p . one can take x(t) having a support so small that 
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a branch of the logarithm log(A J 4(±ir(^ n )e**)) exists near the support of x(i). Then the integrals 
A±(n) are written as 

A±(n) = ^^J e n *^f Alp (±ir(U)e lt ) X (t)dt, 

where, for r/ close to the fixed parameter £, the function &±(j],t) is defined by 

$±(v,t) = log(A A (±ir(7 ? )e it )) - log £>(»/) - irjt. 

We have &±(t], 0) = ±wr/2. The real part of <fr±(r/, i) is non-positive, and is zero if and only if t = 
near the support of x(t). If we put z(t) = ±ir(r/)e**, then the first and the second derivatives are 
given by 



d^ ± (rj,t) = -z(t) 



4 %(.(*)) " 



A^(*(*)) , ^(z(t))A^(t))-A' A (z(t)) 2 



.A A (z(t)) w A A (z(t)) 5 

This shows that, near the support of x(t), dt$±(r], t) = if and only if t = 0. A direct computation 
using (|4.1ip yields 



d ( z ^'a \ _ HW<t>{z) l\a\<P(z) + l 



dz\ X A J~ z(|a|^(z)-l)(|a|0(z) + l) 2 Y \a\<f>(z) - l' 
Substituting z = ±ir(r]) for this identity, we get 



df$ ± ( V ,0) 



2 -|a| 2 



\b\ 

and hence the critical point t = is non-degenerate. Note that each higher derivative of &±(r],t) 
is bounded and \d^±(r],t)\ is bounded from below when (rj,t) is close to (£, 0). This allows us to 
use Theorem 7.7.5 in |H6 to obtain 



±in(n-y n )/2 / m 

A ± (n) = ^ J 1 — — (/^(±ir(6»)) + 0(l/n)) . 

The leading term in the above is a function in £ n . Recall our assumption £ n = £ + 0(l/n). Hence 
the functions in £ n in the leading term, except the exponential term, can be replaced by the values 
at £ of the same functions plus terms of order 0(1/ n). A simple computation shows 

fA*(-ir) = fA*(ir), (r e R\ {0,±1}). (4.13) 
Since e -™(n-y n )/2 _ (•_ 1 )n- 2 /„ e i7r(n- !/ »)/2 ) we get 



e -nH Q (£ n )/2 | 6 | . 

V (i-f 2 )V?--H 3 ( 4 - 14 ) 

x ((l + (-l)«+»»)/^(ir(0) + O(l/n)). 
Next, let us consider the integral Ia2(ji) m f|4. lOj) . Lemma 14.31 shows 



\X A {r{^r l \ < E(V) ■= WA(±r( V ))-'\ = (4 .15) 

v i - V 



for 1/ a/1 + |6| 2 < f] < 1. From this, we have 

\lA2(n)\ < e -"«"l°6r(S™)-logB(C„)) max | 

M=K£«) 
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and hence, by (14.51) . 

e nH Q^«)l 2 \I A2 {n)\ < Ce~ n{losD ^ n) ~ losE ^ n)) . 

It is easy to see that D(rj)/E(rj) > 1 for l/yl + \b\ 2 < rj < 1. Since £ n = £ + 0(l/n), one can find 
a positive constant c such that D(t; n )/ E(£ n ) > 1 + c uniformly in n. This leads us to 

e -nH Q (H n )/2 

I A2 (n) = 0(n 1/2 e- wlog(1+c) ). (4.16) 

\/27rn 

Recalling that ( U n (5o <8> y), ^ ® }e = when n + y is odd, and in view of (|4. 14[) and (|4.16p . we 
get the formula (|4.ip with 

F^O = /^(ir(0)- ( 4 - 17 ) 

4.1.3. Proof of (|4.ip /or i/ie case ^ < in i/ie region (A). In this case, we have — 1 < £ n < 
— 1/ \/l + |6| 2 , —(1 — < r(^ n ) < for every sufficiently large n. We take the circle of radius 

— r (£n) = l r (£n)| for the contour C in (I4.10j) and discuss in a way similar as in the case where 
£ > 0. Note that, when —1 < ij < + |6| 2 , we have 1 7^(77) | = —r[rj) = f{\rj\)~ l . Since 

Aa(.z _1 ) = X A (z), the integral lAi(n) with y = y n \n (|4.10p is expressed as 

By Lemma B~3l we have | AA(?"(|£ n |)e _l *)| < .D(|£ n |). Since -D(r?) = -D(|?7|), the same discussion as in 
the case where £ > leads us to the formula (|4.14p by replacing e J7r ( n ~f™)/ 2 with e J7r ( n -l^«l)/ 2 . For 
the integral lA2(n), remark that Xa (r(v)) = ~~ ^(KM)) for — 1 < rj < + \b\ 2 , and hence, 

by Lemma |4.3[ 

IXAi-rivy^l = lA^KtoDe-*')- 1 ! < IAaCKM))" 1 ! = £(|r/|). 

From this, it is straightforward to obtain 

e nH Q^)l 2 \I A2 ( n )\ < Ce~ n{logD ^ nl) ~ loeE ^ nm , 

which yields the same estimate as (|4.16p . Therefore the formula (I4.ip with the function given 
by (|4.17p is proved. 

4.2. Proof of Proposition 14.11 for the region (B). In this subsection, we assume that the 
parameter £ is in the region (B). 

4.2.1. Some lemmas. Define the function iPb(C) by 

MO = C + i\A-C 2 , CeC\{^R; \x\ >1}, 
where the branch of the square root is chosen to be positive for positive reals. Then tpg is holo- 
morphic in the domain indicated as above, and V'bCC)" 1 = C ~~ iy/l ~ C 2 - The function ipg maps 
the above domain injectively and conformally onto the upper half plane, and it is the inverse of the 
map 4>{z) = (z + z~ l )/2 restricted to the upper half plane. We set, as in (|2.3D . 

X B (z) := ip B (\a\<j)(z)), 

which is holomorphic on the domain Db = D given in (|2.4p . Then the eigenvalues of A(uz) for 
z G Db are Xb{z), Xb(z)~ 1 . Since iPb(—() = — V'-b(C) _1 j we have Xb{— z) = —Xb{z)~ 1 . Note that, 
in this case, IuiXb > 0. For z G Db, define the vectors u B (z),v B (z) G C 2 by the formula (I4.6P with 
replacing Xa with Xb- Similarly, define the functions ipi B (z), ^>2b{z) on Db by the formula (|4.7p . 
Then u B (z) {resp. vb(z)) is an eigenvector of A(ujz) with the eigenvalue Xb(z) (resp. Ab(z)" 1 ), 
and, for any ip G C 2 , the formula (j4.8|) holds by replacing the letter 'A' with 'I?'. Define, for ip G C 2 , 



lAi(n) 



2^ 



-iy n t 



^A(r(\jn\)e 
D{\U) 
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the functions fBip, gBip by (|4.9I) by replacing l A' with l B\ The integral formula corresponding to 
(|4,10p which is obtained by replacing i A' with i B' takes the form 

(U n (5 ®<p),6 y ®?p) £ = oj- y (I B i(n) + I B 2(n)), 
Im(n) = -^ / z-y- l \ B (z) n f B ^(z)dz, I B2 {n) = / ^^-^^(^-"^(z) ds. 

27TZ J C 27TZ J C 

It is easy to see that ub(—z) = —vb(z), v?ib(— z) = —if2B(z), fs-ipi—z) = gs^{z), and hence we 
have 

I B 2(n) = (-l)"+«I B1 (n), (^o®^),^®^)^ = a;-2/(l + (-ir+^J m (n). 

Therefore, to prove Proposition 14. 11 it is enough to consider the asymptotic behavior of the integral 
lBi{ n )- The following lemma is used to specify the solutions of the equation (|4.2p for A = Xb- 

Lemma 4.4. Let rj € M |a| < \w\ < 1. T/ien £/ie solutions of the equation (|4.2p A(z) = 
A_b(z), S -D_b, are given by z = ir(r/) and z = ir(n)~ l , where r(n) is defined in (14, 4p . 

Proof. We have "0s (C) = ~~ iipBiC)/ y/^- ~ C 2 > an d hence 

A' B (z) _ MM*) 



(4.19) 



As(z) " y/l- |a| 2 0(^) 2 ' 

where, as before, /i(z) = (z — z _1 )/2. Then a similar computation as in the proof of Lemma 14,21 
shows the claim. □ 

The following is used to find a suitable path C in the integral formula (|4.18p . 

Lemma 4.5. Let < c < 1 and fi > 0. Let £ = (— 7r < i < ir) be an ellipse given by 

((t) = ccost + i^sini. 

Then we have 

Wb(-w)\ < \MC(t))\ < 

The equality in the first inequality holds if and only if t = — vr/2, and the equality in the second 
holds if and only if t = ir/2. 

Proof. For each (£C\{2;£l; \x\ > 1}, set 

a = Q (C) = Re( v / 1 - C 2 ), P = /3(C) = M v 7 ! - C 2 )- (4.20) 

If C = x + iy, one has a 2 — /3 2 = 1 — x 2 + y 2 , a/3 = — xy. Since a > 0, we see i/jb(() = (a + y)(x/a + i) 
and ^s(C)™ 1 = (a — y)(x/a — i). From this, we have (a 2 — y 2 )(l + x 2 /a 2 ) = 1 which shows that 
a > \y\ and ImipsiC) = a + y > 0. It is straightforward to check 



2a 2 (C) = (1 - x 2 + y 2 ) + y/(l -x 2 + y 2 ) 2 + 4x 2 y 2 , Q = x + iy. (4.21) 

Now, put a(t) = a(((t)). From the facts mentioned above, it follows that |^s(C(*))| 2 = ( a (0 + 
/i sin t) I (a(t) — fi sin t) . Setting 

A(t) :=2W (a(t) _ sin 2 t ) 2 ± mm 2 , 
fi dt 

we see A(t) = 2a(t) 2 cost — 2a(t)a'(t) s'mt. A direct computation using (|4.2ip shows 
cos t 



A(t) = -^P= (1 - c 2 ) y/D$) + (1 - c 2 )(l - c 2 cos 2 t) + ^ 2 (1 + c 2 ) sin 2 t 
where we set 

D(t) = (1 - c 2 cos 2 1 + /i 2 sin 2 i) 2 + 4c 2 // 2 sin 2 1 cos 2 1. 
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By the assumption c < 1, A(t) = if and only if cost = 0. Note that A(t) = is equivalent to say 
that (d/dt)\ip B (((t))\ 2 = 0. On the other hand, we have 

|^b(C(-7t/2))| = = -M + vV + 1 < |Vb(C(tt/2))| = =fi+ vV 2 + 1- 

From this, our claim follows. □ 



4.2.2. Proof of (|4.ip /or £/ie case £ > in the region (B). In this case, we have \a\ < £ < + \b\ 2 

and 1 < r(£) < (1 + |6|)/|a|. We take the circle 

C = C n : z{t) = z n (t) := r(C n )e u (-vr < t < vr) 

of radius r(£ n ) for the contour C in the integrals (|4.18p . Then the curve \a\(p(z(t)) is an ellipse 
since 

\a\(j>(z(t)) = \a\(p(r(£ n )) cost + i\a\(j,(r(£ n )) sin t. 

Since 1 < r(£) < (1 + |6|)/|a|, we have \a\ < \a\(f)(r(£ n )) < 1 and < \a\/j,(r(£ n )). By Lemma [431 
we see 

\^ B (-i\aW{t n )))\ < \X B (z(t))\ = |Vb(|<#M<)))| < |^(i|a|/i(r(C n )))|. 
Note that (±i|a|//(r(^ re ))) = A^(±ir(£ n )). For |a| < 77 < 1, it is easy to check that 



X B {ir{rf)) = W(rj), X B (-ir(r])) = -X B (ir(i])) 1 = iD{n) 1 = i 



l&l-vV 



The integral I B i(n) in (14.180 is written as 

e -nH Q (H n )/2 



lB\{n) 



2vr 



/*7T 

/ e -iynt 


~X B (z n (t))~ 


1 —w 





fBi>{Zn{t))dt. 



(4.22) 



By Lemma [431 when \a\ < rj < + |6| 2 , the equality |A B (r(?7)e rf )| = -D(^) (— tt < i < tt) holds 

en 
h t 

\X B (r(ri)e a )\ 



if and only if t = it/2. Therefore, for a compact neighborhood K C (\a\, + |6| 2 ) of £, there 

exists a constant c > such that 



D(fy) 



< e" c , rjeK, -vr/2 < i < vr/2, |t - vr/2| > e, 



where we fix e with < e < vr/8. Take x(*) G QTW such that = 1 for |i| < e and x(t) = 
for |t| > 2e. By using the above estimate, the integral I B \{n) can be written as 



lBi(n) 
J{n) 



-nH Q (S n )/2 



2ir 



-iy n t 



-(J(n) + 0(e— )), 



A B (r(g n )e^ 



x(t-7r/2)/ B ^(r(&0e a )dt 



-iiryn/2 



X B {ir{in)e l 
D{tn) 



x(t)fB^ir(Zn)e u )dt. 



By (|4.22p . one can choose e > so that a branch of the logarithm log X B (ir(rj)e lt ) is well-defined 
when \t\ < 2e and rj stays in a compact neighborhood K of £. Define 

<&(rj, t) = log \ B (ir(r]) e lt ) - irjt - log D(n), rj E K, \t\ < 2e 

so that the integral J{n) can be written as 

J(n) = e — W 2 /" e «*«».*) x (t)/ B ^(ir(en)e ft )dt. 
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Note that &(r],0) = in/2. By Lemma 14.51 we have Re$(r/, t) < 0. We can also check 

Xb(z) 



■\' B {z) | z \" b {z)\b{z) - \' B {zf 
X B {z) X b (z) 

Thus, by Lemma 14.41 dt&(r/,t) = if and only if t = 0. Furthermore, since, by (|4.19p . 



d ( X' B (z)\ _\' B (z) X" B {z)X B {z)-X' B {zf _ , |a||6| 2 ^) 



dz\X B (z)J X B (z) X B (z) 2 z(l-\a\^( Z y)^ 

we have 

02*, m |a||6|V(r(7/)) (l^fWtEK 

1 {V ' ' (1 + |a| 2 /x(r(r?)) 2 ) 3 / 2 |6| 
Therefore, the critical point t = of &(rj,t) is non-degenerate, and hence Theorem 7.7,5 in |Ho 
gives 



'2vr \b\ 

~(i-e 2 M 

Since £n = £ + 0(l/n), one can replace with £ in the above expression to obtain 



J(n) = a — 1 e in(n-y n )/2 (/ f ir fe )) + 0(l/n)) . 



e -flffQ«n)/2 / |5| 

J B1 (n) = 7 ^^e^-y^\l - 7^%=^ (/B^(ir(0) + 0(l/n)) • 



This shows the formula (|4.ip with 

= /B^(*r(0). (4.23) 

4.2.3. Proof of (|4.ip /or t/ie case £ < m i/ie region (B). In this case, we have —1 < r(£ n ) < 
— (1 — |6|)/|a| for every sufficiently large n. We take the circle of radius |r(£ n )| = r(|£ n |) -1 for the 
contour C in the integral I B \{n) in (|4.18p . Then the argument can be carried out in the same way 
as in the case where £ > 0. Note that X B (\r(^ n )\e lt ) = X B {r{\^ n \)e~ lt ) . By Lemma [43} the integral 
I B i(n) is written as 

p -nH Q (Z„)/2 

Ibi(u) = ( J(n) + 0(e— )) , 



J{n) = / e 



-iy n t 



-it\ 



xit + Tr^f^rdCnir'e^dt 



D(\Zn\) 



where x(t) £ Co°(^) i s chosen in a similar way as in the proof for the case where £ > 0, and the 
function <3?(n, i) is defined by 

$( V ,t) = IogA B (*r(M)e- ft ) - log D(\ V \) - i V t. 

Note that $(n, t) is defined on a neighborhood of (£, 0), and 0) = in/2. By Lemma [4.51 we see 
Re<3?(?7, t) < 0. A direct computation leads us to 

d t $(ri,t) = -i (z^fir+r] 



X B {z) 

i ^'b( z ) i ^"b{z)\ b (z) — X' B (z) 2 . 
d t H V ,t) = -z( J ^ ) +z ^ ) (* = .r(M)e ). 
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From this, it follows that dt&(n, t) = if and only if £ = 0. This critical point t = is non-degenerate 
since 



\b\ 

Thus an application of Theorem 7.7.5 in |Hoj gives 



J <"» = \Zv (1 - a ^-i^ e "° + '" W2 ^ ( - ir(l£l)_1) + ■ 

Noting -r(l^l)- 1 = r(£), we see that the asymptotic formula (|4.ip holds with i^(^) in f|4.23[) 



whenever — 1/y 1 + |6| 2 < £ < — |a|. 

4.3. Proof of Proposition |4"TT1 for the case (C). In this case, we have 



£ = £± := ±1/ v / TTH2, r(£t) = ±r±, r± := (1 ± |6|)/|a|. 

We use the function Ab(z) as in the case of the region (B). But, since z = (1± |6|)/|a| are branched 
points of the function Xb(z), we need to choose a contour C in (|4.18|) so that it does not pass 
through these points. This is done by using Lemma 14.61 below. Let e > be a small constant 
satisfying 



<£<min|l, r+, r?., ^-3 + \J5 + 4r 2 J /2, ^3 - ^5 + 4r 2 J /2| . (4.24) 

Consider the following function: 

z±{t) = z±{r,s,t) := (r| + s) 1/2 cos t + i(r± + r) 1/2 sint, 
where the parameters, r,s,t, satisfy 

{—e < s < for z+(t), , A „ , 

(4.25) 
0<s<e for *_(*). V ; 

Note that if e > satisfies (|4.24|) . we have 



< e < 1, e < r|, y r 2 + e < 1 - e < 1 + e < \jr\ - e, 
and the curve z±(t) is an ellipse in the domain Db- 

Lemma 4.6. Let K C R be a bounded closed interval. Then there exists e > satisfying (|4,24p 

(1) /or every ^ £ K , s £ (—£,0) and \r\ < e satisfying \r — s\ < 2\s\, the function 

log | X B (z+ (r, s, i)) | - £ log \z+ (r, s, t) \ 

has critical points only at t = ±ir/2, and 

(2) for every £ G if, s € (0, e) and \r\ < e satisfying \r — s\ < 2\s\, the function 

log \X B (z-{r,s,t))\ - £log |z_(r, s, t)| 
has critical points only at t = ±ir/2. 
The proof will be given later. 

End of the proof of the formula (14. lj) . Assuming Lemma 14.61 let us prove the asymptotic formula 
(|4.ip for the case £ = £±. First we consider the case where £ = £+. Let {y n } be a sequence of 
integers satisfying £ n := y n /n = £+ + 0(l/n). Let if C (|a|, 1) be a small closed interval centered at 
£+ such that £ n G if for every sufficiently large ra. Take e > as in Lemma 14.61 and fix s G (— e, 0). 



20 



TOSHIKAZU SUNADA AND TATSUYA TATE 



Set r n : = r(£ n ) 2 - r\ . Then r n -> as 7i — > oo. Thus, for every sufficiently large n, the sequence r n 
satisfies \r n — s\ < 2|s|. For any 77 6 .fT, let r^ := r(r/) 2 — r+. Taking .K" small enough if necessary, 
we may assume that \r„ — s\ < 2\s\ and |r„| < e for every r/ £ K. We set 

2(77, i) := z + (r v , s, t) = \J r\ + s cost + ir(rj) shit, — tt < t < tt. 
By Lemma 14.61 the function 

<p v (t) = log \\ B (z(r],t))\ — ?7log |^(77, *)| 

has critical points only at t = ±tt/2. Since 

^(71-/2) = log D{rj) -77 log r(r?) = -H Q {rj)/2 > <P v {-tt/2) = - log D{rj) -rjlogr(r]), 

the function 4> v (t) attains its maximum at t = tt/2 and its minimum at t = —tt/2. As a contour C 
in (|4.18p . we take the ellipse 

C = C n : z n (t) = z(f n ,i) = + s cos t + ir (£ n ) sin t (-tt <t<ir). 
Then we can write the integral lBi{n) as 



2ttz 



By Lemma 14.61 and the discussion above, we have 



z(rj,t) 


-n 




r(r?) 




D(r,) 



V D(in) 



< 1, 



Zn(£) 



(4.26) 



where the equality holds if and only if t = tt/2. Thus one can choose a function x(rj) £ Cq~ 
having a small support such that %(i) = 1 near t = tt/2 and 



Z(77,t) 


-v 


AsOK?7,t)) 


r(rj) 







(??GiT, tesupp (l-x(t))) 



with a constant c > 0. Since Xs(z(r], tt/2)) = iD(rj) and z(t],tt/2) = ir(rj), taking K and the 
support of x(t) small enough, we may assume that a branch of the logarithms log \b(z(t], t)), 
log 2(77, t) exist when t stays in a neighborhood of the support of x(t) an d rj £ K. Then the integral 
isi(n) can be written as 



lBi(n) 



-nH Q (£ n )/2 

2-ni 



3 n*&,t) x( t )f {Zn{t)) ^M dt + ( e - 



where the smooth function $(r/, t), which is defined on a neighborhood of (77, t) = (£,7r/2), is given 
by 

$(77, t) = log As (z (r),t)) - 77 log £(77, i) - log .0(77) + 77logr(77). 

Note that $(77, tt/2) = 7ri(l - 77) /2. By (|4.26p . we have Re $(77,/;) < 0, and the equality holds if and 
only if t = tt/2. The first and the second derivatives of $(77,7;) are given by 



dt$M) 
5 2 $(77,t) 



2 
2; 



2^-77 



,W , {zl)2 y B (z)\ B {z)->: B {zf 



Xb(z) 



Xb{z) 



+ 77|l+(- 
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where we set z = z(jj, t) and z' = dtz(rj, t). By Lemma 14.4} dt&(r], t) = if and only if t = tt/2. By 
a direct computation, we find 



vr/2) 



rl + s 



\a\\b\^(r( V )) 



+ s (1 — r] 2 )^7] 2 — | ct | 2 



r(r?) 2 (1 + |a| 2 /i(r(r/)) 2 ) 3 / 2 r(??) 5 



\b\ 



< 0. 



Note that each higher derivative of $(77, t) is bounded from above uniformly in ij S K and t in the 
support of an d |<9 2< £(r7, tt/2)\ is bounded from below when r] £ K. Since dtz(rj, tt/2)/z(t], it/2) = 

i\/f+ + s/r(r]), Theorem 7.7.5 in [H6J gives 



lBi(n) 



e -nH Q (£ n )/2 2irr(£ n ) 2 



\b\ 



2iri 



n r 



_ e ni(n-y n )/2 % 



r 2 + S 



-fBf(ir(Cn)) + 0(l/n) 



V2 



irn 



\b\ 



\ a-eye-i* 



Mn-y n )/2 {Mir{i+))+0 {l/n)) 



This proves the asymptotic formula (|4.ip with F^{£) = /^(ir(^ + )) in the case £ = £+. When 
£ = take e > and s E (0, e) as in Lemma 14.61 and take the ellipse 



z n(t) = \ f- + s cost + i|j"(£ n )| sint 



for the contour C in (|4.18p . Then, a similar argument as in the case £ = will show the formula 
gU) with = f Bi ,{ir {£_)). □ 



4.3.1. Proof of Lemma \4-6[ Let r, s,t be the parameters as in (14.25p . For simplicity, we write 
z(t) = z±(r, s, t). It is enough to show that the derivative of the function 

f(t) = f(r, s, t) : = log \X B (z(t))\ 2 - flog \z(t)\ 2 , £ G K, 

with respect to t does not vanish when — tt < t < ir, \t\ 7^ tt/2. By a simple computation, we have 



d t f(t) = 2 Re 



Z(t) 



We use the following notations to examine dtf{t). 



We have 



D = D(t) := \z(t)\ = Jr% + s cos 2 1 + r sin 2 t, x(t) + iy(t) := |a|(/>(z(t)), 

fc = (4+s) 1 / 2 , / = (r|+r) 1 / 2 , a(t)+i/3(t) = Vl-|a| 2 ^W) 2 , 
F(t) = (l-x(t) 2 +^) 2 ) 2 +4^) 2 y(t) 2 . 



HMg) * MMg) ., ;2 ,2 

x = — cost, y = — sine, / — k =r — s, 

2a 2 = (1 - x 2 + y 2 ) + >/F, a 2 + /3 2 = VF. 



(4.27) 



(4.28) 
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Then, 

z ^B_ = \ a \ s i n t _ pkfx(D) cos t) - i(ak(i(D) cos t + (3l<p{D) sin t)\ , 

*B D\/F 

■ Sint COSt + I— 7T, 



z 



D 2 D 2 ' 



z' A \ lalfr — s) 
Re ( — J = 1 j^=T sin 1 cos K al< t>{ D ) sin * ~ Pk{Ji(D) cos t) 

+ (akn(D) cos t + Pl(p(D) sint). 

D 3 V i 7 



From this we have 



where we set 



D3 ^ dtf(t) = R(t)(r - s)smtcost + MS(t), (4.29) 
2 a 



for 5(t), we have 



ED\fF 

R(t) = - h al(j>(D) sin t - /3k/j(D) cos i, = akfi(D) cos t + filMD) sin t. (4.30) 

We need to examine the function S(t). Substituting 

xy \a\ 2 kl 

p = = —7r<p\D)ii\D) suit cost 

a aD z 

= k^D)cO S t = ^ _ 2 \^ HD) 2 gin 2 f (431) 

la D z 
Substituting the concrete forms of x(t), y(t) for the second line of (|4.28j) . we have 

I 1 2 / 2 

Alt) =VF + 1- \a\ 2 (j){D) 2 - i^-f- sin 2 t. (4.32) 

Lemma 4.7. There exists e > satisfying (|4.24p such that, for any r,s,t in (|4.25p . one has 

r— |a| 2 Z 2 sin 2 i \a\ 2 l 2 u(D) 2 sin 2 1 1 

> ^2 — + y — + 1 - m 2 ^) 2 - 

Proof. Using the definition of F in (|4,27p . we have 

,2 (/ ^2 W\ 2 l 2 2 V \a\ 2 l 2 fi(D) 2 2 
1- |a| 2 0(,D) 2 + ^ 7i -sin 2 tj +4 1 1 ^ ' sin 2 1, 



from which one has \^F > 1 — \a\ 2 (p(D) 2 + \a\ 2 l 2 sin 2 t/D 2 . To improve this estimate, write 

a| 2 Z 2 2 |a| 2 /V(_ 
^ Sin * + W 



F= l-|a| 2 ^) 2 + ^- 72 -sin 2 t+ 1 1 ^ ' sm 2 t) +G. 



Then, the function G is computed as 

D 2 ~^w — » i ^2 



„ |a| 2 / 2 //(D) 2 ^ . , . 2 2|a| 4 / 2 /i(D) 2 . . . 2 
G = ' ' ^ d(t) sin 2 1 + ' ' 7 ; G 2 (<) sin 2 1, 



where G±, G2 are given by 

L> 2 " zvv L> 2 



d (t) = 2 - |Q|2/ y )2 sin 2 1, G 2 («) = 0(D) 2 - 4 sin 2 t. 
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We would like to show that e > can be chosen so that G > for any r,s,t in (14.25[) . For this 
sake, first note that, since D 2 = k 2 cos 2 t + I 2 sin 2 1, we have 

sin 2 t 111 . . 

For r, s,t satisfying (14.25[) . the following holds. 

fl< (1 + e) 2 < r\ -e < D 2 < r\ + e for z+{t), 

1 < r 2 - e < D 2 < r 2 + e < (1 - e) 2 < 1 for z_(t). 

By this and ()4.33p . G\ and G*2 are estimated as 



G 1 (t)>2-\a\ 2 n[ Jrl±e) , G 2 (t) > 0(D) 2 - 1 > 



Since |a|/i(r-i-) = ±|6|, we can take e > so that 2 — \a\ 2 fi yy r ± e ) > 0> and then G > as 
required. □ 

From ()4.32p and Lemma |4.7[ we see that the function A(t) has the estimate 

i/\ a ,.\ |o| 2 Z 2 w(D) 2 sin 2 i . 
A(t) > A x (t) := 1 1 MV D / + 2(1 - \a\ 2 cj){D) 2 ). 

Lemma 4.8. One can choose e > satisfying ()4.24p w'i/i i/ie property that, there exist positive 
constants A, B such that, for any r,s,t in (|4.25p . we /iave 

Ai(t) > ^l|s|cos 2 t + J Bsin 2 t. 

Proo/. Set f(t) = f(r,s,t) := s cos 2 1 + r sin 2 t so that D 2 = r± + f(t). Since r,s,t satisfy (I4.25p . 
we have \f(t)\ < e < r\. The Taylor expansion gives 

where the constants in the terms 0(f(t) 2 ) can be chosen uniformly in r, s, t. Noting |a| 2 0(?"-|-) 2 = 1, 
we see 

\a\ 2 f(t) / 1\ |a| 2 / 2 . 2 / 1 \ 2 |a| 2 / 2 /(t) . 2 / 1\ r , 
*(t) = (l - ) + sm 2 t(l - ^ + sm 2 * (l " ^ J + 0(/(t) 2 ). 

To rewrite this, note that r_ = r^ 1 and hence 

i-i = ±»!l. i-i = ± 4 ' 6 ' 



9 II ' 4 1199' 

l a\r± r± |a| z r|. 



Substituting this for the above expression of A±(t), we can write 

Ai(t) = l 4 (g±(t) + 0(f(t) 2 )) , g±(t) = TVit) + \b\l 2 S \n 2 t±^P^ S m 2 t. (4.34) 



We need to estimate g±(t) from below. First we consider g+(t). Write 

all 2 . o \ , 9 ( „ ,,9 „ \a\l 2 r 



g+(t) = ( 2 - ^-sin 2 tj \s\ cos 2 t + (J6|Z 2 - 2r + sin 2 t ) sin 2 t. 



r ± 



Then we easily find 
\a\l 2 



sin 2 1 > 2 _ |a|(r + + £) , |6|/ 2 - 2r + MLl sin 2 t > | 6 | r 2 _ (2 _ j 6 | + | a |/ r+ ) e . 



r + r + r + 
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From this, one can choose e > small enough so that there exist positive constants A, B with 

g+(t) = g+(r,s,t) > A\s\ cos 2 1 + B sin 2 t (4.35) 

for any r, s, t satisfying the first line of (|4.25p . The function g~(t) can also be estimated in a similar 
way. Indeed, g~(t) is written as 



g-{t)= ( 2 - sin 2 t \ s cos 2 1 + (\b\l 2 + 2r - sin 2 i ) sin 2 t. 



In this turn, the parameters r, s, t satisfy the second line of (|4.25p . One has 

\a\l 2 r, \a\ Id \a\rl 2 ,., „ / lal la 



sin^ t > 2 - ^ - J-=!-e, |6|^ + 2r - sin^ i > \b\r z _ - 2 + |6| + ^ + ^f e 



Since 2 — jp- > 0, this shows that there exist positive constants A, i? such that (|4.35|> holds by 

replacing g+(t) with g~(t) for any r, s,t satisfying the second line of (|4.25p . Since f(t) 2 = 0{e 2 ), 
the assertion follows from the inequalities (|4.34[) and (|4.35p by replacing g+(t) with g±(t). □ 

End of the proof of Lemma \4-6[ We set 

E(t) = 2aR(t)(r - s) sint + k 2 lfi(D)A(t), 

so that 

2\a\ Jyj 2a 

Since the function a(t)R(t) is bounded for r, s, t satisfying (|4.25j) . we take C > such that 

\2a(t)R(t)\ < C. 

We divide the discussion into two cases. 

(1) Consider the case where s satisfies the first line of (14.25f) . Then we have D > 1 and hence 
//(.D) > 0. Thus there exist positive constants A, B such that 

k 2 ln(D)A(t) > A\s\ cos 2 t + B sin 2 t 

for any r, s, t in (|4.25p . Hence, for r, s with |r — s| < 2|s|, we have 

> (-B- A|s|) sin 2 i — 2C|s|| sint| + A\s\. 

If < e < min{-B/A, AB/(A 2 + C 2 )}, then the polynomial 

p(T) = (B- A\s\)T 2 - 2C\s\T + A\s\ 

is positive for every T, and hence E(t) > for every r, s, t in the first line of (|4.25[) . 

(2) Next, let us consider the case where s satisfies the second line of ()4.25p . In this case, < D < 1 
and n(D) < 0. Thus there exist positive constants A, B such that 

k 2 lfi(D)A(t) < -A\s\ cos 2 t - B sin 2 1 

for any r, s, t. Then, for r, s with \r — s\ < 2\s\, we see 

E(t) < -( J B-^|s|)sin 2 t + 2C|s||sint| - A\s\. 

If < e < min{B/A, AB/(A 2 + C 2 )}, then the polynomial 

p (T) = —(B - A\s\)T 2 + 2C\s\T - A\s\ 

is negative for every T, and hence E(t) < for every r, s,t in the second line of (j4.25p . 
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Therefore, for both cases, taking e > small enough, we see that 

f'(t)= W E(t)cost 

can vanish if and only if cost = 0. This shows Lemma 14.61 □ 

4.4. Functions and Hq(£). We have proved the formula (|4.ip with given by (|4.17|) and 

(|4.23p for the case (A), (B) and (C), respectively. Since 



./fi _ \ n \2 i m 

Xa(«t(0) = As(zr(C)) = iD(0 = 7 =LLp L , | a | < |£| < 1, 

we have u^O"^)) = UB(ir(f)), u A (ir(^)) = v B (ir((,)), and /a^(w*(£)) = /s^ r (£))- Therefore, we 
can write 

*V(£) = ¥>(0(«(0^>. ( 4 -36) 
where the function y>(£) and the vector are given by 

" (f) = r({WCK) + PK)-) ' " (6= ^K)^(f)-'j- ( 37) 

Finally, let us prove the positivity and the convexity of the function Hq(£). 

Lemma 4.9. The function Hq(£) defined by (|1.17p is positive and convex on \a\ < |£| < 1. Moreover 
Hq(£) has no critical points. 

Proof. Since Hq(£) = Hq(—£), it is enough to prove the lemma by assuming \a\ < £ < 1. For 
simplicity, we set <5(£) := Hq(£)/2 so that <5(£) = £logr(£) — logD(£), and hence = logr(£) + 

— %^y- A direct computation shows 

^(0 D'(0 |ft|£ 



-co £>(o (i-ewe-H 2 ' 

Since 1 < r(£) for £ > |a|, we see 5'(£) = logr(£) > 0. It is easy to show 

lim 5(C) = 0, lim $(£) = - log |a|>0, 

£->|a|+0 £->l-0 

and which implies = Hq(£)/2 > 0. We also have 

f (O = -Ro = (i-fVf-l«.|' >0 ' 

which shows the convexity of Hq(£). □ 

4.5. Proof of Theorem 11.41 Taking the modulus square of (|4.1|) with ip = e± and ip = e2, and 
summing up these formulae, we have the asymptotic formula (|1.16p with 

G(0 = \F ei (0\ 2 + l^e 2 (e)| 2 = |^)| 2 (|(u(0, ci )| 2 + |(«(0,ea )| 2 ) = KOI'IKOP, (4.38) 

which is a non-negative function, where the function and the vector u(£) are given in (|4.37p . 
This completes the proof of Theorem 11.41 
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